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I.  INTRODUCTION 


One  aspect  of  designing  an  airplane  or  evaluating  the  performance 
of  a  proposed  design  is  estimating  the  landing  distance  it  will  require. 
In  particular,  some  airplanes  must  land  on  unprepared  fields,  e.g.,wet 
grass  or  ice,  as  well  as  concrete.  The  method  described  in  this  report 
was  developed  to  estimate  lauding  distances  on  these  different  surfaces. 
Sample  results  for  a  STOL  transport  are  presented  and  a  computer  program 
implementing  the  method  is  shown  in  the  appendix. 


2.  ANALYSIS 

2.1  Derivation  of  Equations 
2.1.1  Aircraft  Motion 

The  basic  equation  is  simply 

F  =  ma  (1) 

and  F  is  composed  of  aerodynamic  drag,  the  braking  force,  and  the 
contribution  from  the  thrust  reverser.  Thus, 

F  *  Fdr  +  Fbr  +  Frh  (2) 

where  the  drag  force  is  a  function  of  velocity  and  the  braking  force  is 
the  product  of  the  friction  coefficient  of  the  tires  and  the  downward 
force  on  them.  The  downward  force  in  turn  is  composed  of  the  aircraft 
weight  plus  any  vertical  component  from  the  thrust  reverser.  The  reverse 
tnrust  is  equal  to  a  reverse  thrust  coefficient  times  the  engine  net 
tnrust,  and  net  thrust  depends  on  forward  velocity  and  engine  speed. 

Finally,  engine  speed  varies  with  time.  In  equation  form 


Drag:  Fdr  =  f(V)  (3) 

Braking  Force:  Fbr  -  Cf  *  (GRWT  +  Frv)  (4) 

Horizontal  Component  of  Reverse  Thrust: 

Frh  *  Cv  rev  *  Fgross  (5) 

Vertical  Component  of  Reverse  Thrust: 

Frv  ■  Frh  tan  0  (6) 

Where  0  is  the  discharg  angle. 

Or  >£?'  Thrust:  Fgross  *  f(RPM,  V)  1?) 
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The  forces  are  shown  in  Figure  1  and  the  braking  coefficients  are 
listed  in  Table  1.  Any  moment  generated  by  Frh  is  ignored. 

2.1.2  Engine  Dynamics 

The  engine  is  assumed  to  have  only  two  throttle  positions  - 
idle  and  maximum.  Its  thrust  responds  to  a  change  in  throttle  setting 
according  to  a  first  order  log  equation, 

F  =  Fmax  “  (Fmax  ~  Fidle)  e_t/l  (8) 

where  t  is  the  elapsed  time  in  seconds  and  t  is  a  characteristic  time. 

At  t/r  =  5,  the  transient  is  more  than  99%  complete.  Idle  thrust  is 
assumed  to  be  zero  so  that 

F  -  Fmax  d  ~  e_t/T)  (9) 


TABLE  1 


Braking  Coefficients 


Surface 

Cf 

Ice 

0.05 

Wet  Dirt 

0.15 

Dry  Dirt 

0.25 

Dry  Concrete 

0.30 

2.2  Solution  Method 

2.2.1  Aircraft  Motion 


The  time  step.  At,  is  taken  to  be  0.1  seconds. 


2.2.2  Assumptions  and  Boundary  Conditions 

The  aircraft  is  assumed  to  touch  down  with  an  initial  velocity 
Vo  and  the  engines  idling.  At  time  Tq  (typically  1/2  second),  the  brakes 
are  applied,  the  thrust  reverser  is  deployed,  and  the  throttle  is  opened 
fully.  Full  reverse  thrust  is  attained  after  5  more  seconds,  at  time  T2. 

Some  time  later,  the  airplane's  velocity  drops  below  a  preset  value  V3 
(typically  20  ft/sec)  and  the  throttle  is  closed  and  the  thrust  reverser 
stowed.  This  time  is  denoted  as  T3.  The  brakes  are  applied  until  the 
airplane  is  completely  stopped. 

The  thrust  reverser  applies  both  a  horizontal  and  a  vertical 
force  to  the  airplane,  as  shown  in  Figure  1.  For  the  particular  configuration 
used  as  an  example  here,  the  resultant  reverse  thrust  vector  was  inclined 
at  just  over  57  degrees,  or  very  nearly  one  radian,  from  the  horizontal. 
Consequently,  the  vertical  force  is  not  wasted,  however.  It  increases  the 
load  on  the  tires  and  thereby  contributes  to  the  braking  force. 

Reverse  thrust  loading  is  a  parameter  in  the  analysis.  The 
horizontal  component  is  specified  as  a  fraction  of  aircraft  gross  weight, 
subject  only  to  the  constraint  that  the  resultant  reverse  thrust  not  exceed 
the  total  engine  thrust.  Typical  values  are  in  the  range  0.1  to  0.3. 

Drag  acts  constantly  on  the  moving  airplane.  Its  magnitude  is 
also  a  parameter  in  the  analysis.  The  magnitude  of  the  drag  force  at  the 
Initial  velocity  Vo  is  specified  as  a  fraction  of  the  aircraft  gross  weight, 
rrcm  this,  a  drag  coefficient  is  computed  and  held  constant  for  the  duration 
of  the  integration. 

In  summary  then,  the  landing  scenario  looks  like  this: 

TIME  VELOCITY  EVENT 

0  VO  Airplane  touches  down. 


T1 

VI 

Brakes  applied,  throttle  opened,  thrust 
reverser  deployed. 

T2 

V2 

Full  reverse  thrust  obtained. 

T3 

V3 

Throttle  closed,  reverser  stowed. 

T4 

0 

Stopped . 

2.3  Computer  Program 

ROLLOUT,  the  computer  program  developed  for  this  analysis  is  listed 
in  the  Appendix.  It  is  written  in  small  modules,  each  of  which  does  only 
one  task.  Thus,  it  is  easy  to  change.  For  example,  if  one  is  interested 
in  a  particular  configuration  for  which  he  has  detailed  drag  data,  he  need 
change  only  subroutine  AIRDRAG.  Similarly,  different  engines  may  be  simu¬ 
lated  by  changing  the  tabulated  data  in  subroutine  ENGINE.  ROLLOUT  accepts 
data  in  Namelist  format  as  described  in  the  Appendix.  It  then  computes  a 
set  of  16  solutions,  comprised  of  four  reverse  thrust  loadings  for  each  of 
four  braking  coefficients.  Output  consists  of  a  printed  table  of  distance 
versus  time  for  each  of  the  16  solutions  and  three  different  sets  of  plots, 
two  of  which  are  optional.  The  first  set  of  optional  plots  displays  four 
curves  of  distance  versus  time  corresponding  to  the  four  thrust  loadings  on 
one  plot  for  each  of  the  four  braking  coefficients.  The  other  set  of  op¬ 
tional  plots  displays  four  sets  of  four  small  plots,  each  set  corresponding 
to  a  braking  coefficient.  Each  small  plot  shows  the  relative  contribution 
of  the  brakes,  the  thrust  reverser,  and  aerodynamic  drag  for  the  duration 
of  the  landing  roll.  The  final  plot  is  a  crossplot  of  the  16  solutions, 
showing  rolling  distance  versus  braking  coefficient  for  each  of  the  four 
reverse  thrust  loadings. 

2.4  Results 

The  results  shown  here  represent  a  twin  engine  STOL  transport  weighing 
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160,000  pounds  and  having  a  landing  speed  of  150  feet  per  second  or 
approximately  90  knots.  Braking  coefficients  are  as  shown  in  Table  1  and 
the  reverse  thrust  loadings  were  taken  to  be  0.1,  0.2,  0.3,  and  the 
maximum  possible  (about  0.316).  The  influences  of  the  various  parameters 
are  discussed  below. 

A  set  of  curves  of  rolling  distance  versus  time  (the  first  plot  option) 
is  shown  in  Figures  2,  4,  6,  and  8.  Figures  3,  5,  7,  and  9  show  the 
relative  contributions  of  brakes,  reverser,  and  drag  for  the  same  cases 
(the  second  option).  Figure  10  contains  the  crossplot. 

2.4.1  Braking  Coefficient  and  Reverse  Thrust  Loading 

One  sees  in  Figures  2  through  9  the  perhaps  surprising  result 
that,  except  on  ice,  the  brakes  do  more  to  stop  the  airplane  than  does  the 
thrust  reverser.  This  is  due  in  part  at  least  to  the  particular  reverser 
-sed  here.  The  vertical  component  of  the  reversed  thrust  is  1.557  times 
*be  horizontal.  By  forcing  the  tires  more  firmly  against  the  landing  surface, 
it  contributes  to  the  braking  force.  A  reverser  which  turned  the  exhaust 
jet  more  would  of  course  alter  the  split. 

Figure  10  shows  that  with  the  maximum  reverse  thrust  it  is 
possible  to  stop  the  airplane  on  ice  in  less  than  1000  feet  and  on  dry 
concrete  in  less  than  600  feet.  With  the  minimum  reverse  thrust,  the 
braking  force  is  still  sufficient  to  stop  the  airplane  in  900  feet  on  dry 
cr-’fte,  on  ice  the  distance  required  is  much  larger  than  the  1400  foot 
upper  limit  assumed  (the  actual  distance  was  some  2100  feet).  The  utility 
o:  Figure  10  is  that  is  shows  the  level  of  reverse  thrust  one  must  have  to 
i  .  particular  rolling  distance.  For  instance,  if  the  required 

maximum  distance  is  1000  feet  for  all  surfaces,  the  reverse  thrust  loading 
r.u  t  V  at  lens'  0.3. 
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2.4.2  Engine  Time  Constant 

The  effect  of  engine  acceleration  time  is  almost  impossible 
to  detect  on  plots  of  distance  versus  time.  It  is  a  little  easier  to  see 
in  plots  of  velocity  versus  time,  as  in  Figure  11.  Shown  there  are  two 
velocity  histories,  one  for  a  time  constant  yielding  a  half  second  acceler¬ 
ation.  The  principal  difference  is  a  displacement  of  about  a  quarter  second 
(or  about  9  ft/sec)  and  a  more  pronounced  rounding  of  the  initial  part  of 
Che  curve.  The  sharp  bend  at  the  right  end  of  the  curves  shows  the  use 
of  V3  (here  10  ft/sec)  to  shut  off  the  reversers.  A  cross  plot  of  rolling 
distance  versus  time  constant  is  shown  in  Figure  12. 

2.4.3  Drag 

For  want  of  data  describing  the  aerodynamic  drag  of  the  kind 
of  aircraft  considered  here,  we  used  the  parametric  approach.  For  example, 
DFR  =0.1  means  that  at  the  start  of  the  calculation,  Fdr/GRWT  =  0.1. 
Defining  CD  =  Fdr/V2  gives  CD  -  DFR  *  GRWT/VO2.  Figures  6  through  9  show 
that  for  DFR  <  0.1,  drag  plays  a  very  small  part,  while  Figure  13  shows  what 
can  happen  with  very  large  drag  forces. 
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3.  CONCLUSIONS 


The  computer  program  described  herein  is  a  useful  tool  for  quickly 
obtaining  estimates  of  aircraft  landing  roll  distances  under  various  operating 
conditions.  It  permits  independent  evaluation  of  the  various  parameters 
which  affect  the  rolling  distance  and  can  easily  be  modified  to  model  con¬ 
figurations  other  than  the  twin  engine  STOL  shown. 
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APPENDIX  A 


Input  Instructions 


Rollout  reads  all  its  input  from  Namelist  IN.  All  variables  are 
initialized  in  Subroutine  INPUT  so  that  it  is  possible  to  run  the  program 
without  supplying  any  data. 

The  input  variables  are: 


DFR  The  aerodynamic  drag  force  at  touchdown  (V  =  VO) 

expressed  as  a  decimal  fraction  of  the  airplane's 
gross  weight,  e.g.,  0.1  Default  is  0.0. 


DIS  Logical  variable.  When  set  to  .TRUE,  causes  plots 

of  rolling  distance  versus  time  (the  first  plot 
option)  to  be  created.  Default  is  .FALSE. 


FRAC  Logical  variable.  When  set  to  .TRUE,  causes  plots  of 

relative  contribution  of  drag,  brakes,  and  thrust  re- 
verser  (the  second  plot  option)  to  be  created. 

Default  is  .FALSE. 


GRWT  The  airplane's  gross  weight  at  landing,  pounds. 

Default  is  100,000. 


TAU  Time  constant  for  engine  acceleration,  seconds.  The 

engine  accelerates  from  idle  to  full  power  in  5/TAU 
seconds.  Default  is  1. 


T1  Time  after  touchdown  at  which  brakes  are  applied,  thrust 

reverser  opened,  seconds.  Typically  0.5.  Default  is  0. 

Airplane  velocity  at  touchdown,  feet  per  second.  Default 
is  100.  Setting  VO  negative  stops  the  program. 


VI  Velocity  below  which  thrust  reverser  is  no  longer  used, 

i<t  per  second.  Default  is  20. 
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APPENDIX  B 


Program  Listing 

Program  ROLLOUT  and  all  its  subroutines  are  listed  on  the  following 
pages.  ROLLOUT  uses  the  Calcomp  graphics  subroutine  library,  though  other 
graphics  systems  could  easily  be  substituted.  ROLLOUT  was  written  in 
DCD  Fortran  IV  Extended  and  contains  a  few  dialect  statement  forms  which 
are  not  compatible  with  other  compilers.  All  are  easily  changed,  however. 
There  are  a  few  statements  of  the  form  A=B=C*0.0  and  some  sever 
character  variable  names.  The  forms  READ  *  and  PRINT  *  are  list-directed 
(format-free)  input-output. 

The  listings  contain  occasional  continuation  lines  marked  with  a 
dollar  sign.  These  were  created  by  the  listing  program  to  maintain  the 
right-hand  margin  and  do  not  appear  in  the  actual  code. 


O  r*  ry  r>  n>  oonnnnnn 


F FOG  ° AM  >-Ql.L^«-'T  (TS'PUT  =  i<ii,  OUTPUT,  TA°E 6=0UT°UT  , 
t  CIPT) 
f 

common  /rr\jr^/  VT,  V’,  CPWT ,  THPNAX,  CF,  Tl, 

T  T T Ml,  P VE 
1  (?)  ,  PT,  ?NP,  PFF 

CCHNON  /p.CP/  pr 

PC**HPN  /-.rTA T/  in',  or-F  (250,4),  ORF  (251,0),  fvF 

S  (25P,4) 

PPNMON  /PtrTFl/  rj«,  FrAC 
IfGIPAl  PJS,  FPAC 

PFF  *F;P  THP  BRAKING  PCEFFICIFNtf  FOR  ICE,  WET  ( 

$  pr,y 

CPA^S,  *NP  C FY  PPNCRFTF. 

PTL  TS  TFF  FCFT70NTAI  CCPPONFNT  OF  REVERSE  THRUST 
OTV7ntr  EY  VEHICLE  GFOSS  WFTGHt,  THERE  TS  ALSC  A 
V r t 7 C. l  ccnconfnT,  the  RESULT a NT  PEVFCSF  ■pHFLST 
Vrrior  T?  TNCLINFL  AT  57.  3  CFG  T0  THE  HOFIZON'AL. 

P IHP  NF  T  C  N  CFF  (4),  C  (4,4) ,  NTG  (4),  RTL  (4),  T 

*  (25il )  ,  V  f15 50)  , 

1  Y  (ER'),  X  p  (251,4) 

CATS  pfc  /  <-.f!5,  0.1*,  P.25,  C.?0  / 

PAT A  PTl  /  0.1,  P.2,  0.3,  P.35  / 

PNO  =  n.P 

call  cl p-r  d.r,  i.r,  -  3) 

1  CALL  INFU' 

TF  (VP  .15.  ".()  GO  TO  4 

DFR  TS  THFr  f pA G  F 0°PF  AT  V  =  Vn  AS  A  FRACTION 
OF  GFCSP  WCIGHT . 

OP  JS  ~HE  FOUTVALfNT  CPAG  COEFFICIENT. 

OC  =  CFr.  *  G? NT  /  (V0  *  VP) 

C 

po  x  jrr  =  1,  4 
PF  =  CFE(ICF) 

no  2  TFT  =i,4 

c 

f  RUN  FOUR  ff  VFF  SF  TMFUST  L  OAriNG*  FOP  EACH 

PEAKING  COEFFICIENT. 

c 

CALL  °FR  ^FTl  ( I"’ T ) ) 

PTC  =  F  T  L ( IPT ) 

CALL  rC*‘PLTF  (X,  V,  25P,  NINTO) 

p NO  =  owo  ♦  1. 0 

CALL  CL’TIMJT  (X,  V,  -,  2R",  NINT 6) 


24. 


NTG(I«5T)  r  NIKTG 
0(irF,  IPT)  =  X<  NINTG ) 

no  2  T  =1,  NIVTG 
XF<I,  IFT)  =  XCI) 

2  CCNTINUE 

7F  CAtL  GPAPHIO  (T,  XP,  FJTG|  RTL,  25  lj ,  4) 

IP  (CRAC)  C#IL  QU»D  <PF,  N'T,  FT,  f,  4,  250) 

?  CONTINUE 

C*LL  CROSFIT  (C,  OFF,  ctl  ,  4> 

GO  Tn  i 

4  COLL  SYMBOL  (0.0,  0.5,  0,105,  “FINISHED**,  90.0,  8) 
TALL  FLOTe 
STOP 
FhO 
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c  u  u 


n. 


SUBROUTINE  BRAKES  (COWNr,  tj.,  t,  pf,  FBP  ) 
C 

C  COMPUTES  ERA  KING  FOR^F 

C 

IF  «'r  ,LT.  Tl)  GO  to  i 
FER  =  CF  *  POWNF 

rftupn 

1  FP®  =  P.,1 

MO  BRAKES  FCR  T  <  Ti 


RFTliRN 


TOoonn  ->no 


sifif'nuvTNr  err  (p^t  ) 

/•'HRUST/  FCRC  (E> 

COMPLTrF  ffV.  CO~FF.  FROM  REV.  ThR.  LOADING 

fC*'“r'N  /r7 VFV/  Vf,  V3,  C'RWr  ,  THRMAX,  of,  Tl,  t?, 

S  T  ■*  j  TAU,  CWP 
1  (?)  ,  P^C,  PNC,  DRP 

r  Vc  ( 1 )  =  r""l  *  r,r  WT  /  (1.02  *  f GRO  (1) ) 

CV»<?)  =  1.E577  •  OVP(l) 

1  .5R77  =•  T^(?7,’  OFT.  ) 

1.0?  =  ?*C.«** 

GUApp  «C-«INC1  RfPULTftNT  TV  *  1. 

T  F  (rvrM1)  •  L c  •  ".5UP3)  PfTl'PN 

p 

r  thppp  ARE  T WF  CPP^NE  ANC  CINE  Dr  57.3  OEG . 

C 

rVR(l)  =  E<*C3 

rvF(?)  =  C.PM1P. 

pTL  =  *  FGRTd)  f  GPWP 

r. 

C  i.*37<*  =  l.R?»**.543J 

C 

PFIV)  *,  ~rTL  FAD  7"'  RE  ffCITEP  To  rtl 

PPTU  ^N 

PNn 


n  ^  o  ti  o  o  ^non 


RtlBFOUTlN*  CCi'FUTF  (X,  V,  T,  NO,  NINTG) 

FCMMON  /GIVEN'/  VC,  V2,  G*WT  ,  THPVAX,  CF,  Ti,  T2, 

?  H,  Tflu,  TVP 
1  (?) ,  PTC,  PNC,  OFP 

COMMON  /FETAPC/  TUT,  RCF  <?50,4>,  ORF  (25", 4),  RVF 
«  <?*0,4) 

TNTFGPATES  THE  FOLIATION  OF  NflTIFN,  F  =  MA  ,  IN  A 
TAYLOP  SEPIFS  FXPANSION,  VIZ 

X (T+GT )  =  X (T)  ♦  OTMCX/OT)  ♦  ( CT**2/2) *( C2X/CT 2) 
WHFPF  DX/HT  =  v 

02x/riz  =  *  =  f/n 

F  =  TFAG  ♦  PRAfTNG  PFV.  THRUST 

“CGIN 


r 

C 

r 


OTMFNSTON  T  (NO),  V  (NO),  X  (NO) 

rT  =  o.i 

oT2  =  cy  *  m  /  2,r 

xm  =o.o 

V(l)  =  vr 

T (1 )  =  0.C 

OH  =  32.174  /  GPWT 

T?  =  36C3.0 

T-»  =  3FC%0 

NM»y  =  NC  -  2 

BPF(1,  TPT)  =  0.0 
PFF(1,  IPT)  =  0.0 
RVF(1,  IR T )  =  0.0 

TF  (OFF  ,Gt.  c#n)  p’F(l,  TFT)  =  l.n 

tntfgpaticn  LOOP 

DO  3  I  =2,  N**X 

CALL  AlFTPAf  ( V ( T  -  1>,  TP»G) 

FALL  FNGINF  (V(T  -  1),  ',>o*AX,  r t MCRA5) 

CALL  SFCOLUF  (T(l  -  i),  Tl,  T2,  T3,  TAU,  THRMA>, 
*  THP) 

CALL  RFITHP  (ThP,  CVF,  V3,  V(T  -  1),  FPH,  FR  V) 
FALL  “FAKES  (GPWr  ♦  FP  v  »  Tl,  t(j  _  *)  ,  CF,  FBP) 
TF  (FPH  .GT.  O.n)  FRH  =  FPH  ♦  favdfAG 
PSUK  =  (“RAG  ♦  FPP  ♦  FPH 

TF  (FROM  .GT,  0.  0)  GO  TC  1 


B,p ( T , 

TFT) 

S 

BPFd  , 

TP“r) 

ORF  (T, 

TFT ) 

s 

nPF ( 1 , 

IPT) 

F  V  F  (  I  , 

TRT) 

T 

P  VF ( 1 , 

TPT) 

GO  T0 

2 

BP  F  (T  , 

TFT) 

r 

FPP  / 

FRUM 

OR  F ( I , 

IPT) 

s 

OF  AG  / 

FSt'M 
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o 


c 

c 

c 


c 


c 

e 

c 


NTN1G  TS  TK’  NUMTR  OF  INTGRflrTON  STE°S. 

U  NTNTf,  =  T 

IF  <?3  .GF.  ’er-i.r)  »3  =  T(T) 

rfTurN 

p  *<n 
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ooo  ooo  r>  n  .t  oo  n 


subroutine  crnspLT  (hist,  cff,  ftl,  n) 

CCHMON  /FRA Mf /  XF,  YF 

DIMENSION  5C  (6),  CFF  (N),  OIST  (N,N),  RTL  (N)  ,  X 
S  (6) 

DAT*  YNUM  /  6.0  f 

plots  landing  piptan,'e  vs  prakilg  coefficient 

WTTH  PFVEPGF  r  HP JST  LOADING  «S  A  PARA  MeTEP • 

DATA  X  (5),  X  <  6»  /  4C0.0,  200.0  / 

PATA  V*'  (5),  6C  (6)  /  0.0,  / 

XF  =  A.O 

YF  =  5.T 

PALL  PI  0T  (2.0,  2.  n,  -3) 


CALL 

axis  O 

.  t- ,  0.0,  13HPRAKING  COFF., 

-  11» 

XF, 

$  O.n, 

ac  ( 5 ) , 

1  PC 

(6)  ) 

CALL 

AXIS  (0 

•  T,  O.n,  8HPIFTA  NCE ,  6,  YF, 

9  3  •  ^  « 

X  (  E) 

3  X(F>)) 

csll 

cyhsol 

(3.5,  L.q,  '.  i ,  thftl,  ij.o, 

3) 

YC 

= 

YNUM 

PO 

1  T 

=  1,  A 

PC  (I) 

=  CFF(T) 

1  CO 

NTINUP 

03 

3  I 

=  1,  N 

no  ?  j  =  i,  u 

X (j)  =  Pisr ( J,  I) 

2  CONTINUE 

LA LL  Ct IF  (EC,  X,  N,  1  ,  0,  O 

Y  LCCAT T CK  T  CR’M  PEVpPSr  THR'JFT  LOAOING  VALLE 

YN  =  (X(A)  -  x*  ■*  0  •  C  >  /  2"f'.r. 

MA<t  ? Uk E  THEY  D3N*T  OVpplAF 

IP  <<YO  -  YM  .IT.  0.3)  YN  =  VO  -  0.3 
YO  =  YN 

G  *  LL  NU'EER  (3.5,  YN,  0.1,  RTL(I),  0.0,  3) 

3  CONTINUE 

CALL  PL  O'  (/>,(  ,  -2.  f» ,  •?) 

I  FT  Ufc‘ N 

f  *L,r 
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O  O  C-* 


c 

c 

r 

r 

C 

C 


* 

« 


SU!SDOU'‘TNc  rFG  TWF  <V,  Thr*-AX,  *AVopAG) 

^CRMON  /'‘FPUST  /  FGSL  (  E  ) 

OIMFNGTO*’  Rm  (5),  V<?L  (6) 

nrr*  MCI  /  t,r.,  311. 64,  167.46,  223.28  / 

Dfi‘rf  FGSL  /  42780.3,  4CPPC  .%  403PC.O,  5C  000.r, 

El’QfiP.o  / 

PATA  Knct  /  rtCt  Z60C.r,  5inn.fi,  770J.O,  187QC.C 
/ 


r'ATA  F''t=  OK*-  r  NG7  F'E  *  T  ERA  LFVFL. 
cGSl  GROSS  th?UST  (Cv  =  1.8). 
rnsi.  j s  tHC  p;k  DRAG. 

V  tk  vft.  -rA  °  l  E 


1 

2 


IF 

(  V 

.  GT  .  VSl (3)  ) 

GO 

TO 

2 

IF 

(  V 

.OT.  VSl  (2) ) 

GO 

TO 

1 

I 

=  1 

GO 

T0 

4 

I 

=  2 

GO 

TO 

4 

TF 

(V 

.GT.  VSl (4) ) 

GO 

TP 

y 

▼ 

r  7 

GO 

TO 

4 

I 

=  4 

«  PSU*<r  -  |>r  c  KG  TME r 


4  7 

THrMfl  X 
*  PGEL  (TH) 
TAKHRAG 
t  ’?ncL  (I)  )) 
PPTUPN 
FNn 


(V  -  VSL  <  I>  >  /  (VSMT  4-  1)  -  VSL  ( I )  > 
2.“  *  (FGSUI)  ♦  7  *  (FGSLCI  ♦  1)  - 

2.*  •  (PDSl  CT3  4  I  •  IPDSLtl  4  1)- 
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l>  o  o  o  o  <->  o  -)  o  o  o  ->  ■>  r>  r>  r> 


SUP:<OUTJN'"  GSAFhtf  (T,  X,  N f  G>  r-r  L  ,  No,  MC) 

COMMON  /GIVFN/  V  fi,  \/3,  GPW  ,  rHP“AX,  OF,  Tl,  T2, 
?  1  3,  i flU,  CVR 
1  (2),  KTO,  RNC,  rFr 

COMMON  /FkAMF/  XF,  YF 

OTMENSI0K  N*o  (VO),  r-'TL  (MO),  T  (NO),  X  (NO, MO) 
OAT*  YNUM  /  5.  7  / 

ON  (  INC  CALCCMf 

D°flW  AXFS  AND  PLOT  X  VS,  T 

XF  =  5.r 

YF  =  6.T 

CALL.  PLOT  (C.C-,  1.5,  -3) 

CALL  AXIS  (P.'S  P.n,  4HTIMF,  -  4,  XF,  J .  0 ,  0.1, 
f  4,0) 

CALL  AXIS  (P.r,  0.",  aHHISTANCF,  8,  YF,  90.0,  0. C, 

♦  200.1) 

A  NNul A  “F 


CALL 

SYMBOL 

(  0  •  3r  , 

5.8, 

0.  14, 

1"2, 

0.0, 

-1) 

CALL 

SYMBOL 

(0.45, 

5.6, 

0.10, 

1H=, 

0.0, 

1) 

CALL 

NUM6FP 

(0.65, 

5.8, 

J.1P  , 

CF,  0 

.0, 

2) 

CALL 

S  Y  V  oOL 

<0.30, 

5.6, 

0.10, 

5HCF* 

~  , 

0.  0, 

CALL 

NU*5FF 

( r.9o , 

5.6, 

u.  IP, 

OFF, 

c.o, 

2) 

nvAw  eox 


CALL 

PL  0" 

(0. 2f, 

5.5, 

3) 

CALL 

PLOT 

(0. 20  , 

6. 11 , 

2) 

CALL 

PLOT 

(1.35, 

6.0, 

2) 

C*LL 

PI 

(1.35, 

5.5, 

2) 

oAl  L 

PL  CT 

( b • ?P  , 

5.5, 

2) 

CALL 

SYMBOL  (4.5, 

5.  9 

r 

9 

YC 

=  YNUM 

00 

l  T 

=  It 

MO 

Nt.M 

=  NTG(I) 

♦  l 

NS 

=  NM 

♦  1 

3HPTL, 


0.0, 


3) 


I  NS*  PT  SCALE  F  A  CT  OPS  INTO  T  flPc-SY. 

TN  !  =  T  (NH) 

INS  =  T  (NS) 

T(NM)  =  t1.0 

1  (NS)  a  4.0 
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x(w*,  i)  =  ?or.r 

r.M.L  CLIc  X(1,T>,  ktg<  r  > ,  1,  o,  p> 

v  ( =  *  m- 

’(VS)  =  NS 

V,J  iS  HE  V~Ft  ICA  L  LOCATION  OF  '  Hr  =>'|_  W9tl 

YN  =  X(N7G(J),  I)  /  ?C0.<' 

jr  ( (yo  -  YN)  ,Lt,  0,3)  YN  =  YO  -  0,3 
YO  =  YN 

TILL  Nl'Herp  ((4,5,  yn,  0.1,  FTL(T),  0.«,  3) 

CONTINUE 

C^LL  DL  0T  (ft.*,  -1,5,  -X) 

kE'  U"N 


SUSFOUTTNF  ikfut 

PCMHON  /GIVEN/  VP,  V 3,  GFWT  ,  th9max,  CF,  Tl,  T2, 

$  T3,  tau,  rvR 
1  (?)  ,  RTC,  PNC,  OFF 

TCMHON  /PLCTFl/  PIS,  F^AC 
l^GICAl  PIS,  FP'C 

HAT  A  VO,  V  3,  GFWT,  TmFNAX,  PF,  "1,  T3,  TAU,  CVP  / 
«  1P0.0,  2*\C, 

1  i.OEOF,  S.^F  C  <*,  C.  2,  '’.E,  O.C,  1.0,  Cl.  5,  0.1  / 
data  Tfr  /  r.r  / 

DATA  OTS,  FF  Ac  /  2  *  .FALSE.  / 

NAME! 1ST  /  IN  /  DFP,  CTS,  F^AC,  G^WT,  TAU,  Tl,  VO, 

J  V  3 

PFAP  IN 

RETURN 

FND 
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?UpF'OUT  TNp  mpUT  (X,  V,  7,  NO,  NINTG) 

CCWMON  /rivr*v  VP,  V 3,  GFVT,  -wp*'AX,  OF,  Ti,  T2, 

?  7  3,  TfiU,  ovn 

1  (?) ,  prp,  RNC,  PFp 

c 

r  O^INT  Hr f OFF 

C 

0TMpNSipK  -  (NO),  V  (NO),  X  (NO) 

P*TNT  Gr  WT,  THRMAX,  CVR,  RTP,  OFR,  Tl,  T  2,  72, 

*  V3»  VO, 

1  nr,  -|A.u,  RNC 
C 

C  PFTNT  V,  X  ARRAYS 

r 


no  2  K  =  3,  3 
PRJNT  lnl 

Ti  =  1  +  10  C  *  (X  -  1) 

T?  =  T]  f  49 

T?  =  )TNP(T?,  no  -  ?,  NINTG) 

nO  1  T  =  II,  72 

°rtnt  in?,  t  a )  ,  v  (i) ,  x  (i) 

TF  (I  +  5P  .IF.  NINTG)  ppjnt  in'*,!  d  ♦  5P),V  (I 

3  ♦  50) , X  (7  + 

1  5") 

i  C0NTT*'Uc 

TF  (7?  ,  rf,  MNTG  ,  OR,  Tl  +  lPP  ,FT,  NINTG)  P=T'JFN 
Pp  I NT  ire, 

?  CONTINtiR 
RfTUF  N 

12  0  FOREST  (  1H1 ,  7X,  *GR.NT.  **  F8.P,  *  MAX  THRUST 

♦  **  r»,c, 

1  *  OV7W  =*  rr.-»,  *  CVF V  =♦,  F5.3,  *  RTL  =  *  F5.3, 

*  ♦  pFP  =*  F  5, 3  f  AX,  *’l  s* 

c  Ff> .  2,  *  T  2  =*  F6«  ? ,  *  T?  =  ♦  F5»  2 ,  ♦  V3  =*  F7,(, 

•*  *  \jr  =• 

3  p7.7,  *  OF  =*  F 5, 3 ,  «  TflU  =*  F5.3,  *  RUN  NO.  s* 

*  F’.n  ) 

101  fopmjt  (  i HO ,  ? (6X,  AH'r INF ,  2X,  8HVELOCITY,  ?X, 

♦  8  HPI P 7  ANr  F ,  lTX) 
i  ) 


F3P*  *  T 

C  1H  , 

3F1P. 

?  ) 

r %r 

(  1H», 

ACX, 

3FJ  n,  ?  ) 

FC^H AT 

<  i  H 1 

) 

cun 
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o  n  n  n  o  .">00000 


SURPOUTINF  OUAC  (rc,  NTG,  P~L,  T,  N,  NT) 

COMMON  /RITAFC/  7?T,  ^PF  (250,4),  DPF  (250,4),  R \JF 
i  (250,4) 

COMMON  /PRAME/  XF,  YF 

riMj£N<ION  NTG  (N),  RTL  (N)  ,  T  (NT),  XO  (4),  YO  (4) 

PLOTS  f  GFAFH  FO*  PACH  OF  FC'JP  "  HRUST  LOADINGS. 

-  A  CH  G"APH  SHOWS  THF  PELATIVE  CONTRIBUTION  OF 
S  OPAG, 

BRAKES,  ANC  THFUST  pPVFRSFR. 

OATA  XO  /  2.0,  4.0,  -4.0,  4.0  / 

P  AT  A  YO  /  5.0,  0.0,  -5.0,  C.O  / 

DATA  FI,  F?  /  0.75,  1.0  / 

YF  =  4.f 

XF  =  3.0 

CALL  PL  O’-  (l.C,  1.5,  -3) 

CALL  FACTOR  (FI) 
nO  1  L  =  1,  4 
C*LL  PLOT  (XC(L),  YOU),  -’> 

K'L  =  NTG  (L ) 

OO  2  T  =  1  ,  NL 

RVF(X,  L)  =  R  VF ( I  ,  L)  +  9RF(I,  L) 

ORP (I  ,  L)  =  DKF ( 7  ,  ))  ♦  r VF ( I ,  L) 

2  CONTINUE 

NM  =  NTG(L)  ♦  1 

NS  =  NM  +  1 

'  NM  =  T  (NM) 

TNS  =  T(RP> 

T(NM)  =  O.fi 

T  (  NS)  =  5." 

RRF(WM,  L)  =  CF  F  ( N  M  ,  L)  =  ®VF(N«,  L)  =  0.0 
BRF (NS,  L>  =  OPF(NS,  l)  =  <VF(NF,  L)  =  0,25 

i.kLL  axts  ( r . a,  r.o,  4htim-,  -c,  xf,  3.0,  0.3, 

1  5.0) 

~AcL  AXTS  ('’.'>,  fiHFFACTTON,  6,  YF,  90.0,  0.0, 

?  *.25) 


CALL 

Cl  IP 

a, 

3  R  F  ( 1 , 

L>  , 

N’ G(L3  , 

1, 

o, 

0) 

C6  LL 

Cl  TP 

<7, 

h  VF(1, 

D  , 

NT  G (L )  , 

1, 

o, 

0) 

CALL 

CLTP 

<T, 

C  K  F  (  1 , 

L), 

N'r  G(L) , 

1, 

r, 

0) 

T  (  NM)  =  TNM 

I  (NS)  =  TNS 

annotation? 

1HRUST  LOFTING 

CALL  SYMBOL  (1.9,  0.C5,  0.1",  5H»TL  =,  3.3,  5) 

CALI  NUMEFP  (2.5,  0.C5,  0.10,  PTL(L),  C.O,  3) 
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<"»  o  O 


ViJL  =»  =  2.0  *  BRF  ( 3?  ,  L)  -  0.05 

Y'LS  =  ?.T  *  (8-F(3P,  L)  ♦  ’VF(38,  L))  -  C.T5 

C*-LL  ^Y-'BCL  (0.5T,  YBL6 ,  0.10,  BHBkAKES,  0.0,  F) 

^  a »  L  SYMBOL  (0.5P,  YTl.B ,  P.ir,  7HTHP  RFv,  0.0,  7) 
COLL  FY  BOL  (0,15,  1,  PL,  0.  IP,  4HDPA&,  3  0.0,  4) 

1  r04TXMUf 

BkAKlNG  COFFFIrIENT 

call  sr-'eoi  (-i.pc,  -0.75,  p.20,  102,  0.0,  -i> 

FUL  SYHaOL  (-P.75,  -0.75,  O.lf,  1H=,  P,(,  1) 

CALL  NU“F“P  (-r.cp,  -p.75,  P.1C,  OF,  C.O,  2) 

C-LL  F*  C” or  (F2) 

CALL  PLC'r  (5.C,  -1.5,  -3) 

KF7  LT’N 
FK'P 


o  o  n  r>  r> 


SU8P0UTTNP  PFVTHP  rrHP,  CVR,  V3,  V ,  FPH,  FRV) 

rOHPUTrs  HCFT7CK,’AL  ANO  VrPTTCfL  COh=>ONENTS 
nr  PrvFPSE  thrust. 

OTMrNSIPN  rvp  (2) 

IF  (V  ,LE.  V3)  GO  TO  1 
FPH  =  the  *  CVPti) 

FFV  =  THR  *  rvF (2) 

RFTU°N 

HO  PFVFCSE  THPUST  MHFH  V  <  V3 

1  FPH  =  FPV  =  o 

PFTUPN 
FNn 
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ono  r> 


SUP*  OITTNF  FFfoi  UP  (T,  '1,  T?,  T3,  TAU,  TH?MAX, 
f  THP> 

C 

C  rOMFL'TPS  FKGINF  "HRUST  pr’-M£CW  JDLF  AND  PULL 

C  USTK'G  *  FTPST  CPnc-R  i_sg. 

C 

IF  ("  .IT.  ti  ,CF.  T  ,GT.  T?)  GO  "0  l 

"HO  =  tffwx  »  (1.0  -  PXP<  -  (T  -  Tl)  /  TAU) 

r 

FULL  POVF«  IF  «?F»CHn  AT  T  =  T?. 

IF  (<"  -  ’1)  /  ~*U  ,GF.  5,r  .ANP.  T2  .GE.  3600. C) 

?  T  2  =  " 

PFT'JPN 

T  HROrT  l  F  T5  rLOSFC  F?«?  T  <  7 1  OF  T  ►  T  3 
1  THO  r  f,r 

PFTUFN 
Fwn 
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SUBROUT  IK 

P 

Cl  ?F  (X,  Y,  N,  K 

COMMON  /7.r  ZYZN/  rl  .  S2,  Lc 

LOGICAL 

Si,  s? 

otme-nstok 

X  (1),  Y  (j) 

XTNtXI 

r 

(X  -  XMTN)  /  rx 

YIN (  YT 

= 

(Y  -  YMTN*  /  CY 

XNTN 

r 

X«N  ♦  11 

YMJN 

r 

Y(N  ♦  1) 

OX 

r 

X(N  ♦  2) 

OY 

= 

Y(N  ♦  2) 

XI 

= 

XTMX(l)  ) 

Yl 

= 

YIK(Y  (1)  ) 

si 

= 

J  .GT.  r 

S? 

= 

•  FPLSF, 

LS 

r 

L 

1  T 

=  2,  N 

X2 

=  XTN(XfT)) 

Y2 

=  YIN ( Y ( I  )> 

TF  (J  . 

GT  • 

*)  S2  =  BOO  (  T , 

r  A  LL  AGAINST  (XI,  Yl,  X2,  Y2) 
XI  =  X2 

Yl  =  Y  2 

SI  =  .FALSF. 

c2  =  .FALSE. 

1  CONTTNUr 
PFtupn 
FMO 


,  U 


EO«  0 
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Sli^°OUTTKc  (XI,  VI,  X2,  V?) 


r 

c 

c 

f 


r 

c 

r 


f 

r 


CLTpc  ^hc  y  ikp  p: qw  (XI, Vi)  T0 

.1  <=  x  <=  xr,  r  <=  v  <=  yf 


r  CMM^m 

/FFftMr/  xp, 

Yc 

FfMHON 

/  7S  ZY  ? v  /  cj 

.  S ?, 

L CGI  CM. 

si,  <•?, 

TS 

X? 

=  xCCF(X2, 

Y2> 

X  i 

=  xrrp  (X< 

,  YI) 

IP  (XI 

,P0.  0  .avn. 

K2  .FC. 

T )  GO 

TP  (  (XI 

•  SNT,  X  ? )  . N  F •  *>) 

pr 

T  (JF'N 

Je  (  K  1 

,a,r.  *n  gp  Tri 

SWAP  fkjp$ 


T 

= 

XI 

Yl 

= 

X? 

X? 

r 

T 

T 

= 

Yl 

Yl 

r 

Y2 

V? 

r 

▼ 

K 

r 

XT 

r 

x? 

X? 

r 

X 

T  C 

r 

FI 

<1 

r 

F? 

^2 

s 

■*■  c 

TF 

(  (  <1 

•  C 

•  J.)  . 

Lcpt 

r  t  p  £ 

(X  = 

*> 

r 

C 

r 


r 


r 


Yl 

=  Yi 

- 

XI  * 

(Y2  - 

Yl) 

XT 

-  c. 

r 

GO 

-n  i 

IF 

(  (XI  . 

A  KG  • 

2) 

.  FQ. 

P)  GC 

T  A 

RIGHT 

FTCF 

(X 

=  XF) 

Yl 

=  Yl 

♦ 

(XF  - 

Yl  )  « 

>  (Y2 

XI 

=  XF 

nr 

T  0  1 

IF 

(  (  XI  . 

nn, 

t) 

.  PO. 

(t)  r.  c 

TC  5 

P(1  PK 

(Y  - 

T) 

1 

1 

=  XJ 

- 

Yl  * 

(Y2  - 

XI) 

Yl 

=  n. 

»> 

4? 


(X2,Y2)  AGAINST 


*0  6 


(X2  -  XI) 


Yl)  /  (X’  -  XI) 


(V?  .  Yi) 


n  n  o  n  o 


go  To  1 

5  IF  ((<1  ./INC.  £>  .  F  Q.  P>  GC  T0  1 

r 

c  "">p  (Y  =  YF) 

c 

XI  =  XI  ♦  (YF  -  YD  •  (X2  -  XI)  /  ( Y2  -  Yl) 

Y 1  =  YF 

G  r  T  o  l 

IF  w€  PFACh  HFFE,  Twr  | TNT  F*?ON  (XI, Yl)  TO  (X2, 
*  Y2) 

TS  VISIPLF. 

6  CALL  Pt  CT  (XI,  Yl  ,  31 
CALL  PLCT  ( X 2 ,  Y2,  2) 

IF  (SI)  CALL  SYHBOL  (XI,  Yl,  0.105,  LS,  C.C,  -1) 

TF  (S2)  CALL  SYMBOL  (X2,  Y2,  0.105,  LS,  0.0,  -1) 

F FT  URN 
FNO 


